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Abstract 



We develop Truncated Conformal Space (TCS) technique for perturbations of c = 1 Conformal Field The- 
^-j- . ories. We use it to give the first numerical evidence of the validity of the non- linear integral equation (NLIE) 
derived from light-cone lattice regularization at intermediate scales. A controversy on the quantization of Bethe 
states is solved by this numerical comparison and by using the locality principle at the ultraviolet fixed point. 
It turns out that the correct quantization for pure hole states is the one with half-integer quantum numbers 
originally proposed by Fioravanti et al. Once the correct rule is imposed, the agreement between TCS and 
NLIE for pure hole states turns out to be impressive. 



p; 1 Introduction 

cu : 

The scaling functions of 1+1 dimensional integrable models on the cylinder have proven to be 
a very useful non-perturbative tool of investigation of their finite size effects, in particular of 
the renormalization flow properties of the vacuum as well as of the excited states. Various 
methods have been proposed for the calculation of these important quantities, as the Truncated 
Conformal Space (TCS) method [||] or the Thermodynamic Bethe Ansatz (TBA) 0,0]. A very 
promising method is the one introduced some years ago by Destri and de Vega |4], |] (similar 
methods were independently introduced in Condensed Matter Physics by other authors ||). 
It consists in defining a lattice model and give evidence that its continuum limit reproduces 
the sine-Gordon QFT (sG) 0, then deriving a non-linear integral equation (NLIE) which is 
basically the continuum limit of the Bethe equations. There is one such equation for each 
excitation in the spectrum of the physical states. The problem is to give a general rule to write 
down the NLIE corresponding to a state characterized by a certain distribution of Bethe roots 
and holes. This problem, namely to find the NLIE for excited states, was first addressed, in 
a QFT context, in [§] which discussed the case of pure hole excitations. The general setup 
of excited states was described in detail later in |J. The solution of the NLIE, the so-called 
counting function, is a central object in Bethe Ansatz approach to integrable QFT from which 
it is possible to reconstruct the eigenvalues of all the local conserved currents of the theory put 
on a cylinder of circumference L. 
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The aim of this paper is to present a numerical comparison between the NLIE and TCS. 
We develop a TCS for c = 1 CFT perturbed by its cos (3(j> operator, that defines the sG model 
as perturbation of a CFT. We compare TCS data against scaling functions computed from 
numerical integration of NLIE and find a very good agreement, especially in the attractive 
regime, but also in the repulsive one. This agreement, however, only shows up for a specific 
choice of the quantization rule for the Bethe Ansatz states which is different from the one 
reported in || and instead agrees with the one made in ref.||. This leads to the observation 
that the NLIE produces in fact more states than those present in the sG Hilbert space. It is 
known that the relation among the Hilbert space of c = 1 CFT, the sG model and the massive 
Thirring (mTh) model is a quite delicate issue |l(J. It is out of the scope of the present letter 
to investigate this very important question. We intend to return to it in a more extensive paper 
where a more detailed analysis of the results presented here will be given PH . It is important to 
keep in mind that the NLIE (and in general Bethe Ansatz methods) are only able, at present, to 
reproduce states with even topological charge (i.e. even number of solitons minus antisolitons) 
in sG model. 



2 The NLIE and its properties 

The Lagrangian of sG theory is 



C = Jdx Q<9^$ + ^ : cos -)j 



Often it is convenient to use the parameter p = or the 6- vertex anisotropy 7 = Q 

We do not reproduce here the deduction of the NLIE from the lattice Bethe equations. 
Although there are subtle remarks to do on all the derivations presented so far || ||, this aspect 
is out of the scope of this short letter and will be dealt with in a forthcoming longer publication 
[Tl| . Here we are content to present the final equation and explain how it can be used to 
investigate the finite size corrections to energy-momentum of the excited states. 
The NLIE for the sG/mTh model reads 

/+00 
dxG{\ -x-irj) log(l + (_i)<V^+^)) (1) 
-oo 

where 77 is a positive free parameter, small enough not to reach any singularity of the counting 
function Z off the real axis. With such prescription, eq. ([]]) is totally independent of the choice 
of r) by Cauchy's theorem. I is related to the circumference of the space-time cylinder L (the 
spatial volume) by I = ML, where M is the mass scale of the theory, chosen to be equal to the 
soliton mass. The kernel G(X) is given by 

/•+ 00 sinh (p ~ 1)fc7r 

G(\)= / dke* kX 2 



sinh E*E cosh *f 



1 In terms of p, p = 1 is the free fermion point, and p = 4 , A; = 1, 2, . . . are the thresholds where a new 
breather appears, p < 1 corresponds to the attractive and p > 1 to the repulsive regime. 
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The source term g(X\Xj) is specific to the state considered and vanishes in the case of the 
vacuum, which in Bethe Ansatz language corresponds to a sea of real roots. 
The counting function Z(X) has the property that any Bethe root satisfies 

e iZ(x 3 ) = ( 2 ) 

The converse is not true, i.e. there are values A^ for which e %zlyXh) = — (— I) 5 that are not 
solutions of the Bethe equations. Such locations on the real axis are called holes. 

Excited states are seen, in a typical Bethe Ansatz spirit, as different root configurations for 
a system of Bethe Ansatz equations and are given by creating some holes in the vacuum sea 
and adding some complex roots. Holes and complex roots are 0(1) in number, while the real 
roots of the vacuum sea are 0(N) (N being the number of sites), and tend to infinity in the 
continuum limit. The number 5 is equal to the number of roots modulo 2 on the lattice. 

For a state having H holes at positions A^ on the real axis, N c complex roots A c (of dose type) 
such that l^mAcI < 7rmin(l,p), N w complex roots (of wide type) with |3mA c | > 7rmin(l,p) 
and N s special roots/holes A s : 

H N c N w N s 

g(M^j) = E - - E ^ A - Ac ) - E *"( A _ ^w) _ 2 E *( A - A *) 

h=l c=l w=l s=l 

where Xj denotes the set of all parameters (j = h, c, w, s). The function x(A) is the odd primitive 
of G(X) and 

Xu(ty = xW + sign(p — 1) x(X — in min(l,p) sign(Q : m A)) . 

is the so-called second determination of the function X- The interested reader can find more 
details in 0. 

The parameters Xj are determined by the Bethe conditions 

Z(X 3 ) = 2nI J , I je Z + ^ (3) 

which are the continuum analogs of the logarithms of the Bethe equations. Any root of the 
Bethe equations, as well as any hole, satisfies (||) and its position is therefore fixed by giving its 
Bethe quantum number Ij which is half-integer for 5 = and integer for 5 = 1. 

The number of special roots/holes X s , for which Z'(X S ) < 0, as introduced in |§, need 
not to be specified a priori for a state. Indeed one can easily convince oneself that in the deep 
infrared limit / — > oo, Z\X) is always positive and these objects never appear. Their appearance 
in some states for values of / smaller than a certain l is dictated by the fact that sometimes the 
logarithm in the convolution term of the equation fails to remain in its fundamental domain, 
and an analytic continuation imposes to modify the equation itself by adding new terms to the 
source. The number H e ff = H — 2A^ S is invariant for a given state as the spatial volume / varies 
and coincides with the number of dressed solitonic particles in the state. 

The number of the various species of roots and holes in a configuration is not arbitrary but 
rather is subject to constraints originating from an analysis of the asymptotics of the lattice 
counting function. These constraints can be recast in a form that does not contain the number 
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of real roots, and therefore their validity can be extended to the continuum limit, in a form that 
we call the counting equation |J 

H eff = 2S + M C + 2M W 0(p - 1) 

where 0(x) denotes the usual step function and S is the total spin of the Bethe Ansatz system, 
which can be identified with half of the sG topological charge Q. 

Once the function Z(\) is evaluated as a solution of the NLIE, it can be used to compute the 
energy-momentum (E, P) (or in principle the eigenvalue of any other local integral of motion) 

/+oo J\ 
—e ±x log(l + (-1)V Z(A+ ^) 
-oo 27r 

where 

/ N h N c N s 

z(X\Xj) =M ij2 e±Xh ~ e±Ac + Q ( 1 ~P)^2 e±Xw i 1 + e Ti ^ +1 > isa ^ mX ^) - 2 e±Xs 

\h=l c=l w s=l 

In particular, P ± (L) for L — > oo reproduces the asymptotic particle content of the state. 
For L — > instead, it is expected to be related to the central charge c and conformal dimensions 
A± of the ultraviolet CFT state by P ± = -^(c - 24A ± ). 

Notice that eq. (HD does not coincide with the one reported in ||. The argument of one 
of logarithmic terms in || differs by a phase (— l) 5 . Our choice gives the correct infrared 
asymptotic 0(e~ l ) for the IR asymptotics of finite size effects, as it must be in a massive 
theory. The expression in would instead give a rapidly oscillating integrand, which gives a 
contribution that can be (after tedious calculations) proven to be only 0(1/1) in the IR limit. 
A careful analysis of logarithmic branches confirms that our choice is the correct one. We shall 
report the details about this point in a later publication [TT|. 



3 UV behaviour of multi-soliton states 

We consider sG theory as the perturbation of a c = 1 massless free boson $ compactified on 
circle of radius R by a potential V = ~ : (e l/3 * + e~ % ^) :, where g = The interaction 

term has conformal dimensions Ay = Ay = @- = and becomes marginal when (3 2 = 8tt 
which corresponds to p = oo. With the convention^ adopted here, the primary states of the 
unperturbed CFT, have left and right conformal dimensions 

A , 1 f n 1 \ 2 n 2 m 2 R 2 nm 

A* = - — ± -mR = — - H ± 

n ' m 2 \R 2 J 2R 2 8 2 

They are created by the vertex operators V^ m , where n e Rand the winding number m (which 
is the eigenvalue of the topological charge Q) satisfies m G Z. Locality of the operator product 

2 Note that this convention is different from the one usually adopted in the literature (see |jl^| and references 
therein), where the compactification radius is 1/R and the elementary sine-Gordon field is the dual of the 
conformal scalar at the UV fixed point. In our convention, the conformal scalar field is the UV limit of the 
elementary sine-Gordon field. 
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algebra puts a severe constraint on the operator content of the UV field theory. At a generic 
value of the coupling constant (3 there are only two different possible maximal local subalgebras 
of the algebra generated by the vertex operators ||10), namely a bosonic algebra A b generated 
by operators with {n, m G Z} which corresponds to the modular invariant partition function 
and is the UV limit of the sG model, and a fermionic algebra Af generated by operators with 
{n G Z, m G 2Z} [J{n G Z+|, m G 2Z + 1} corresponding to the T 2 invariant partition function 
and the UV limit of mTh field theory. The two algebras overlap in the sector with even winding 
number (topological charge) which is the sector common in the two off-critical field theories and 
is the one that so far has been studied using the NLIE. We limit ourselves to this sector in this 
paper too. 

The ultraviolet behaviour of the solutions of the NLIE can be obtained in the so-called a kink 
limit" well-known from TBA. The method of calculation is the same as outlined in [|]|f§ so ^ 
us just give the results. For the vacuum (quantized with half-integer Bethe quantum numbers) 
one gets the leading behaviour 



with c = 1, which is the correct UV central charge. 

For a two-hole state quantized with integers (which is the prescription of 0) and having a 
left-moving hole with quantum number I + > and a right-moving one with quantum number 
I~ < (with 1+ ^ I~) we obtain 

A ± = ^- + -R 2 --±I ± , 
8R 2 2 2 

which corresponds to the vertex operators V±± 2 and their descendants. As these operators are 
not included in any of the two possible maximal local operator algebras, we have to exclude the 
integer-quantized two-hole states. The UV dimensions of the same two-hole configuration with 
half-integer quantization turns out to be 

A ± = -R 2 - - ± J± 
2 2 

which instead describes Vo,2 and its descendants, the primary state corresponding to the minimal 
choice ^ = This is in fact the vertex operator that one expects to be the UV limit of 
the lowest lying two-soliton state in view of the results in |10|. The results presented in || 
are different when specialised to the two-hole states. Since (as shown in the next section) 
the numerical results confirm our calculations we think that the calculation of the ultraviolet 
dimensions in [|] should be taken with care. 

For generic states with even number of solitons if, half of which move to the left and half 
to the right, with their quantum numbers not necessarily symmetrical to the origin, choosing 
half-integer quantization we obtained that the UV limit corresponds to the family Vq } h (the 
primary state given by the lowest possible choice of Bethe quantum numbers). This contradicts 
the rule in || which is integer quantization for H/2 odd and half-integer for H/2 even, but is 
in accordance with what is expected from the UV conformal field theory and agrees with the 
numerical results we present in the next section. Integer quantization, instead, leads to the 
vertex operators V ± i H and their descendants, which are again excluded by locality. 
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4 Truncated Conformal Space at c = 1 



In this section we describe results obtained using the Truncated Conformal Space (TCS) method 
which support the validity of the NLIE describing the vacuum and pure hole states. 

The TCS method was originally created to describe perturbations of Virasoro minimal mod- 
els in finite spatial volume flTJ. We have developed an extension of the method to study per- 
turbations of a c = 1 compactified boson, more closely the perturbation corresponding to sG 
theory. 

The perturbation keeps the winding number m and the momentum P conserved. Therefore, 
the Hilbert space of the massive theory can be split into sectors labelled by the values of P and 
m, which are quantized by integers. The TCS method consists of retaining only those states in 
such a sector for which the eigenvalue of the UV conformal Hamiltonian He ft is less than a 
certain upper value E cut , so the truncated space is defined as 



H TC s(s, m, E cut ) = : P |*) = s , Q = m |*) , H CFT |*> < E cut |tf )} . 

For a given value of s, m and E cut this space is always finite dimensional. In this space, the 
Hamiltonian is represented by a finite size matrix whose explicit form is: 

H = T {L„ + L --I + 9 —^Bj , 

where L and L are diagonal matrices with their diagonal elements being the left and right 
conformal weights, I is the identity matrix and the matrix elements of B are 



B 9>v = ($| : cos(/3y?(l, 1)) : |$) . 
We choose our energy units in terms of the soliton mass M, related to the coupling constant g 



by the mass gap formula given in ref. 12 



\ 2-2Av 

The TCS method provides a nonperturbative way to obtain the spectrum (the mass gap, the 
mass ratios and the scattering amplitudes) numerically which can serve as a tool to check the 
predicted exact results of integrable field theories and get a picture of the physical behaviour 
even for the non-integrable case. The systematic error introduced by the truncation procedure 
is called the truncation error. It increases with the volume L and can be made smaller by 
increasing the truncation level (at the price of increasing the size of the matrices, which is 
bound from above by machine memory and computation time). 



7lT(l - Ay) I 



6 



While we do not describe technical details of the TCS computation (we refer the interested 
reader to QTJ| ) , we do have to make some remarks on how the TCS method applies to c = 1 
theories. 

The (dimensionless) energy of any state \1/ (the scaling function) goes with the volume L as 



E*(L) 7r(c-12(A+ + A^)) 



M 
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Bl + J2Ck{*)l k{2 - 2Av) 



(4) 



k=l 



where A^ are the left/right conformal dimensions of the state in the ultraviolet limit, B is the 
universal bulk energy constant (the vacuum energy density) and the infinite sum represents the 
perturbative contributions from the potential V. 



The bulk energy constant has been predicted from TBA and reads [ 1 2 



B = -~tan(^) (5) 



4 V 2 

(the same was obtained from the NLIE in ||). This is a highly non-analytic function of p and 
it becomes infinite at the points where p is an odd integer. In fact, at these points there is a 
value of k for which k(2 — 2Ay) = 2, and Ck ( 1 i r ) — > oo. The infinite parts of B and Cu (\&) 
exactly cancel, leaving a logarithmic (proportional to I log/) and a finite linear contribution to 
the energy, by a sort of a resonance mechanism. 

In the attractive regime the results of TCS converge very well when increasing the truncation 
level E cut (as far as we do not approach too much the free fermion point p — 1). On the contrary, 
in the repulsive regime there are convergence problems for the scaling functions E^(L). These 
are related to certain ultraviolet divergences in conformal perturbation theory already known in 
the literature (for a detailed discussion see |13|), which disappear when we subtract the vacuum 
energy from the excited states energies (which is a sort of a renormalization procedure). The 
relative energies £^{L) = E^(L) — E vac (L) converge well, which can be seen e.g. from the fact 
that the mass gap acquires its exact value (which is 1 in our units) quite rapidly when increasing 
the truncation level. Consequently, in the repulsive regime one can only trust the relative scaling 
functions produced by the TCS method, while in the attractive regime we will see that even 
the absolute energy values agree very well with the predictions of the NLIE (including the 
predicted bulk energy constant (H) , which is completely analytic for p < 1 and thus logarithmic 
corrections are absent as well). 

Even in the attractive regime, to get reliable results we had to work with dimensions around 
4000, which was just around the limitation of our machines. Diagonalization for one sector 
took around one day for twenty evenly spaced values of the volume L. This shows that the 
TCS for c = 1 theories is far less convergent than the one for minimal models (in the original 
Yang-Lee example the authors of took a 17 dimensional Hilbert space (!) and arrived to 
very accurate results). Such a slow convergence meant that the total procedure (building up 
the basis, calculating the matrix elements and diagonalizing the Hamiltonian) had to be done 
completely in a compiled C program. For the matrix diagonalization, we used the LAPACK 
numerical linear algebra library, freely available from Netlib over the Internet. 
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Figure 1: The first few energy levels in the vacuum (m = 0) sector at p = | (plotted with 
diamonds) for -E cut = 17.0 (the dimension of the space is 4141) and the NLIE prediction for the 
vacuum scaling function (shown with a continuous line). 

5 Numerical results 

5.1 TCS and NLIE in the attractive regime 

We have made the comparison of the numerical data from TCS and the NLIE predictions at 
several values of the parameter p. For illustration, we present the case of p = |. For this and all 
other values of the coupling constant, we obtained a spectacular agreement between the results 
obtained by the two methods, up to deviations of order 10~ 4 — 10~ 3 which are smaller than the 
size of the points in the plots. The deviation grows with the volume L, exactly as expected 
for truncation errors. By studying other parameters such as the mass gap, the breather-soliton 
mass ratios and the rate of convergence of the energy levels with increasing the value of E cut , 
the small differences can be clearly attributed to the inaccuracy in the TCS data. 

Figures 1, 2 and 3 show the results in the m = 0, 2 and 4 sectors, respectively. The value of 
the predicted bulk energy (|]) has been subtracted from TCS data in order to normalize them 
the same way as the NLIE data are. Note that in order to achieve agreement the holes should 
be quantized with half-integers. We have made some trials of integer quantization in the NLIE, 
but the results obtained differed significantly from the TCS data, with errors of the order 10%. 

5.2 The repulsive regime 

In the repulsive regime the absolute scaling functions do not converge well in the TCS method as 
it was described in Section 4 above. We take as an example the value of the couplings = 1.5 and 
plot the relative scaling functions obtained from TCS, comparing them with the corresponding 
ones obtained from the NLIE (figures 4 and 5). One can see that the truncation errors become 
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Figure 2: The first few energy levels in the m = 2 sector at p — | (plotted with diamonds) 
for _E CMt = 20.0 (the dimension of the space is 3917) and the NLIE prediction for the two hole 
scaling functions with quantum numbers (— |, |), (— |, |) and (— |, |) (shown with continuous 
lines). 
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Figure 3: The first few energy levels in the m = 4 sector at p = | (plotted with diamonds) 
for E cut = 26.0 (the dimension of the space is 3403) and the NLIE prediction for the four 
hole scaling functions with quantum numbers (— §, — |, |, |), (— |, — |, |, |) and (- 
(shown with continuous lines). 
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Figure 4: The first few energy levels relative to the vacuum in the m = 2 sector at p = 1.5 
(plotted with diamonds) for E cut = 20.0 (the dimension of the space is 4445) and the NLIE 
prediction for the relative two hole scaling functions with quantum numbers (— §, |) and (— §, §) 
(shown with continuous lines). 

larger, at values of the volume I close to 10 the deviation can be observed even from the figures. 
For I < 5 the agreement is still within an error of order 10~ 3 . Once again, results obtained with 
integer quantization for the holes showed a large deviation from the TCS data. 

6 Conclusions 

This letter should be seen as a preliminary report of the work done by our group to better 
investigate and clarify some of the properties of the NLIE in sG theory. We have given the first 
numerical check against Truncated Conformal Space data - the only checks so far have been 
calculations of the conformal spectrum at small scale (kink limit) and comparison with the 
Factorized Scattering Theory at large scale. The results presented here are the first evidence of 
the validity of NLIE in the intermediate range of scales and thus show that the NLIE describes 
correctly the interpolation between the UV and IR asymptotic regions. Our approach has also 
resolved some ambiguities about the choice of quantization rule for the Bethe quantum numbers: 
pure hole states are always quantized by half-integers in the sG model on a cylinder with periodic 
boundary conditions. Work is in progress to extend the results to complex roots, the outcome of 
which will be published - together with details omitted in this paper - in a forthcoming article 

E • 
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Figure 5: The first few energy levels relative to the vacuum in the m = 4 sector at p = 1.5 
(plotted with diamonds) for E cut = 22.5 (the dimension of the space is 4149) and the NLIE 
prediction for the relative two hole scaling functions with quantum numbers (— |, — §)§,§), 
(— §, — f j f) |) an d ( — §> — |j |) |) (shown with continuous lines). 
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